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ABSTRACT 

The remarkable similarity between emission spectra of narrow line regions (NLR) in 
Seyfert Galaxies has long presented a mystery. In photoionization models, this similarity 
implies that the ionization parameter is nearly always the same, about U ~ 0.01. Here 
we present dusty, radiation-pressure dominated photoionization models that can provide 
natural physical insight into this problem. In these models, dust and the radiation 
pressure acting on it provide the controlling factor in moderating the density, excitation 
and surface brightness of photoionized NLR structures. Additionally, photoelectric 
heating by the dust is important in determining the temperature structure of the models. 
These models can also explain the coexistence of the low-, intermediate- and coronal 
ionization zones within a single self-consistent physical structure. The radiation pressure 
acting on dust may also be capable of driving the fast (~ 3000 km s -1 ) outflows such as 
are seen in the HST observations of NGC 1068. 

Subject headings: galaxies: active — galaxies: Seyfert — ISM: general — line: formation 



1. Introduction 

Optical spectroscopy of the ionized gas around galactic nuclei provides strong constraints on 
the excitation mechanisms. Line ratio diagrams of active galactic nuclei (AGN), such as those 
of Veilleux &: Osterbrock (1987), allow us to distinguish three major classes; nuclei excited by 
starbursts, and the two types excited by a bona-fide active nucleus: the Seyfert narrow line regions 
(NLRs) and extended narrow line regions (ENLRs), and the Low Excitation Nuclear Emission Line 
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Regions (LINERs). As far as the Seyfert galaxies are concerned, the compilations of Veilleux & 
Osterbrock (1987), Veilleux(1991a,b,c), and Veron-Cetty & Veron (2000) show that the emission 
line ratios are remarkably uniform within a given object, and very similar (with typically less than 
0.5 dex variation) from one object to another. 

The standard paradigm proposes that the NLR and the ENLR are excited by photons originat- 
ing at or near a compact nuclear source (see, e.g. Osterbrock (1989)) having a smooth featureless 
power-law, or broken power-law EUV ionizing spectrum. Within this model, the clustering of the 
observed line ratios within such a restricted domain of parameter space presents a problem, since in 
such models the excitation is in large measure controlled by the ionization parameter, U, and this 
is free to vary over a wide range. The implications of this is that there should exist photoionized 
regions with very high excitation corresponding to high U, but objects such as these are missing in 
the Veilleux & Osterbrock (1987) diagnostic diagrams. Instead, the constraints derived from the 
line ratio diagnostic diagrams imply that U has a similar value in all NLRs. Modelers have been 
therefore forced to make the arbitrary (and possibly unphysical) assumption that the gas density 
in the ionized clouds must fall exactly as the inverse square of the distance from the nucleus. 

One rather successful variation on the photoionization models is to assume that the power- 
law EUV spectrum passes through a highly-ionized screen of diffuse gas, before being absorbed in 
the vicinity of dense (much higher pressure) clouds (Binette et al. (1996), Binette et al. (1997)). 
Although such models produce generally good results when compared with observed ratios, these 
have the difficulty that the geometry of the ionized gas is ill-defined and necessarily becomes a 
new free parameter. These models can also be questioned on physical grounds, since the assumed 
pressure difference between the cloud and inter-cloud components is not explained. 

Another variation is to consider the integrated emission from a whole family of photoionized 
clouds of different densities (Komossa k, Schulz (1997); Ferguson et al. (1997)). These models can 
be made to fit the observed spectra quite well, but the geometry of the clouds and their relative 
contributions to the emission remain free parameters and are poorly constrained. 

In any of these photoionized models, the (often quite violent) internal dynamics of the ENLR 
are not addressed. These include clear evidence of non-gravitational motion, and evidence for 
outflow at velocities up to (or even in excess of 1000 km s" 1 ) {eg Pedlar et al. (1989), Allen et al. 

(1999) ). It was results such as this that led Wilson & Willis (1980) to suggest that the nucleus 
ejects radio components that interact with ambient gas and replenish the high kinetic energy and 
ionization of the NLR. This led to the development of photoionization by fast radiative shocks, which 
generate a strong internal photoionizing radiation field (Dopita & Sutherland (1995), (1996)). This 
model is successful in reproducing the strong lines observed in many AGN (see review by Dopita 

(2000) ). In addition, it can explain other features which simple photoionization models fail to 
address. For example, where the NLRs are spatially resolved, there often are found to be strong 
correlations between radio power and either line luminosity de Bruyn & Wilson (1987) or line 
width Wilson & Willis (1980). Such correlations exist not only for Seyfert galaxies, but persist up 
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to much more luminous classes of radio galaxies. However, the fast shock model has its limitations. 
It has difficulty reproducing the very strong coronal lines of [Fevn], [Fex], [Svill], [Sxn], [Sivi] 
or [Sivil] (Wilson & Raymond (1999)) which imply photoionization with very high local ionization 
parameters. These models also cannot reproduce the large outflow velocities inferred in the [O in] 
line in some objects such as NGC 1068 (e.g. Kraemer & Crenshaw (2000); Cecil et al. (2001)). In 
shock models, most of the [O m] emission is expected to be generated in the photoionized precursor 
of the shock, which theoretically should be the undisturbed gas in the galaxy. Furthermore, the 
NLR lines of (probably) the majority of Seyfert galaxies are quite narrow and consistent with line 
broadening due to galactic rotation rather than being produced in shocked outflows. 

In this paper we argue that some of these properties can be explained in terms of radiation 
pressure dominated dusty photo-ablating clouds. In the past, dust has been considered in photoion- 
ization models both through its contribution via photoelectric heating (e.g. Bottorff et al. (1998)) 
and as shielding component thanks to its large opacity in the FUV. However, the effects of dust 
photoelectric heating and of the radiation pressure on dust on both the ionization and hydrostatic 
structure of the NLR clouds has not been examined in a systematic fashion before. In the next 
section we detail these effects and show how these provide a natural explanation as to why the 
inferred NLR ionization parameter seems to be confined to the range —2.5 < log U < —2. In §3 we 
describe the basics of our dusty, radiation pressure dominated models. These are compared with 
other types of photoionization models, allowing us to separately understand the effects of photo- 
electric heating and of radiation pressure on the emergent spectrum. The predictions of all these 
models are compared with observations of NLR of Seyferts on some key emission line diagnostic 
diagrams. The concluding remarks are given in §4. In following papers we will investigate further 
the effects of different parameters on our radiation pressure dominated photoionization models and 
present a full hydro dynamical model exploring such structures. 



2. A New Paradigm: Dusty, Radiation-pressure Dominated Photoionization 

The effects of radiation pressure from line absorption have previously been considered in pho- 
toionization models of the emission line regions of AGN (Elitzur Sz Ferland (1986)). This is im- 
portant in determining the internal density structure at large ionization parameter. The effect of 
dust within the ENLR has also been considered previously (Netzer & Laor (1993)). However, the 
radiation pressure exerted by the ionizing source acting upon dust within the photoionized gas has 
usually not been accounted for by the standard photoionization models, with a few notable excep- 
tions (e.g. Baldwin et al. (1991)). The standard photoionization modelling assumes an isochoric 
(or constant density) structure within the ionized slab instead of an isobaric (constant pressure) 
structure. Dust is important not only because of the radiation pressure acting upon it, but also 
because of the photoelectric heating which it produces in the photoionized plasma, and in the way 
it competes with the gas in absorbing the ionizing photons. 

However, all of this critically depends on whether dust can survive the either the strong EUV 
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radiation field, or the relatively energetic conditions characterizing the photoionized region. In 
the ionized plasma of the ENLR, grains are electrically charged, locking the motion of the dust 
to that of the ionized gas. Both the ionized gas and the dust is likely to be streaming at roughly 
sonic velocities from ionization fronts located around photo-evaporating dense clouds located in the 
"ionization cones" of the AGN. The problem of dust survival is therefore reduced to the question 
of whether it can survive for a time equal to its passage through the high-emissivity region near 
the cloud surface as discussed in §3. Our detailed photoionization models show that the average 
gas temperature is ~ 10 4 K, at which the thermal sputtering of grains is negligible and incapable of 
destroying the grains in so short a timescale. This is in contrast with the case of radiative shocks 
in which gas temperatures reach ~ 10 5 — 10 6 K and grain destruction by sputtering occurs on very 
short timescales (Tielens, McKee, Seab, & Hollenbach (1994)). Similarly, our models show that 
grain temperatures are too low for thermal evaporation of grains by quantum heating provided 
that the dust grains are large enough (> 3 nm ). A third potential dust destruction mechanism is 
that of Coulomb explosion, in which the grain charging becomes large enough that the Coulomb 
repulsion between the atoms composing the grain exceeds the atomic bond strength. According to 
our models this condition is met for grains smaller than about 10 nm when the ionization parameter 
is of order unity. In the models described here we have simply set the minimum grain radius equal 
to 10 nm to allow for the destruction of the smaller grains. We will return to the detailed physics 
of grain destruction processes in these environments in a future paper. 

In NLR regions observed by Oliva et al. (1994), Moorwood et al. (1996), and Kraemer & 
Crenshaw (2000) the coronal lines appear remarkably strong as compared to H/3. The existence of 
such coronal gas requires that both the ionization parameter and the radiation pressure are very 
large (Binette et al. (1997)). Elements such as Fe and Si are normally largely sequestered in the 
grains. The strength of the coronal lines of these elements requires that in the coronal gas, at 
least, the dust has been largely destroyed. Our model, as depicted in Figure 1 is therefore one 
of a dusty, radiation-pressure dominated region surrounding a photoevaporating molecular cloud, 
which in turn is surrounded by a coronal halo within which the dust has been largely destroyed. In 
this paper, we will consider only the emission spectrum generated by the dusty dense inner region 
near the stagnation point in the flow around the cloud. This will be the highest emissivity region 
in the flow, and it will be the region in which the radiation pressure gradient is matched by the gas 
pressure gradient, allowing us to construct static rather than dynamic photoionization models. 

Consider So> the ratio of the available pressure of ionizing radiation, $, to the gas pressure at 
the irradiated photoionized surface of a photoevaporating cloud 1 : 



where k is Boltzmann's constant, c the speed of light, n and T e are the hydrogen density and the 
electron temperature at the irradiated cloud surface (indicated by a subscript zero elsewhere). We 



^Prad $o/c 



~ Uq 




(1) 



Pgas 2.3nkT e 



1 This definition implies a factor 2.3 lower value than that of Krolik et al. (1981) 
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define 5* = $0/ (eeuv) = / to be the ionizing photon flux (cm _2 s _1 ) impinging at the 

surface such that Uq = S*/nc is the canonical ionization parameter. {e EVV ) (= {hu EVV )) is the local 
average EUV photon energy. 

In hydrostatic equilibrium, the gas pressure gradient must match the local radiative volume 
force exerted by photon absorption, and this force can be separated into its two components, the 
force acting on the gas through photoionization, and the component due to dust absorption; 



where tp{v) is the local radiation flux at frequency u, n^X^ 1 ) the number density of ion +i of 
atomic species m, a v {X+?) the corresponding photoionization cross section with threshold v l m . In 
the dust absorption term, « a b s {y) and K sca (y) is the dust absorption cross section and scattering 
cross section respectively, normalized to the hydrogen density and g(v) = (cos0 u ) is the mean 
scattering angle of the dust grains at this frequency. The boundary condition is given by Pq, the 
pressure at the irradiated face. 

We can readily estimate a critical value of Uq above which dust absorption dominates over 
the photoelectric absorption. For an equilibrium photoionized plasma, the absorption of ionizing 
photons in a slab is simply equal to the local recombination rate in the plasma: 



where is the recombination rate to excited states of hydrogen. The absorption of ionizing 
photons by dust is given by: 



where k is the effective dust opacity of order ~ 1CP 21 cm 2 per H atom. Dust absorption becomes 
relatively more important as the strength of the ionizing field increases, and dust becomes the 
dominant opacity in the plasma when: 



For a typical value of as = 2 x 10 13 cm 3 s 1 , the critical ionization parameter is ~ 0.007. 

This value is closely related to that at which radiation pressure starts to dominate either the 
gas pressure or the dynamical acceleration of the plasma. If (e E uv) is the mean energy of the 
ionizing photons absorbed by gas, and (e FU v) is the mean energy of the ionizing absorbed by dust, 
equation (2) can be approximated by: 





dS^/dx = -asfi 2 



dS*/dx = —nnS* 




(3) 



UBn 2 (e E uv) + -S* (e FU v) nn 



■ (£euv) £ TT 

c ip 



(4) 



-6- 



where (e EVV ) / (e FVV ) = ip , and the effective ionization parameter for the dust-absorbed photons 
U FVV = £U. Generally speaking, both £ and ip will be greater than unity. Clearly, dust becomes 
the dominant term in this equation when: 

TpQB 

Uo > -z — 5 

£ck 

this expression is the same as equation (3) except for the ratio ip/£ , which is of order unity. In the 
typical ionization field expected around an AGN, (fo(v) oc v~ x ' and with a 'standard' dust model 
with solar abundances, a Mathis, Rumpl & Nordsieck (1977) (MRN) size distribution n(a) cx a -3 ' 5 
and a lower size cutoff of 10 nm , we estimate that this condition gives a critical ionization parameter 
U ~ 0.01. 

If the ionization parameter into a static plane-parallel slab of ionized gas is increased above this 
critical value, then the radiation pressure gradient imposed as the EUV radiation field is absorbed 
through the ionized slab will induce a corresponding local gas pressure gradient. Thus, if the initial 
(outer layer) pressure of the ionized gas is Po and the radiation pressure in the EUV field is P ra d 
then the gas pressure close to the ionization front is Pq + P ra d- Thus, at high ionization param- 
eter, the pressure in the ionized gas close to the ionization front is determined by the externally 
imposed ionization parameter, Uo, and the local ionization parameter becomes independent of the 
initial ionization parameter. In this way, the emission line spectrum in the low- and intermediate- 
ionization species will become effectively independent of the external ionization parameter. This is 
the core of our explanation of why the emission line spectra of NLR are so similar one to another. 

While Uo defines the ionization parameter at the face of the cloud, one can define a mean 
ionization parameter U = S*/nc = Uono/n where n is simply the mean density, that is an emissivity 
averaged density using a weight oc n 2 e . Binette et al. (1997) showed that U tends asymptotically 
towards a constant value in the limit of high Uo- 



3. Results of Models 

We have constructed static, plane-parallel, radiation-pressure dominated, dusty photoioniza- 
tion models using the code MAPPINGS III, which includes the implementation of eqn. (2) to 
evaluate the local radiation pressure gradient and hence the local density gradient. MAPPINGS 
III treats several of the properties of dust such as photoelectric heating, as described previously 
by Dopita & Sutherland (2000). The dust model includes both silicaceous and carbonaceous grains 
with a MRN grain size distribution and uses the absorption and scattering data of Laor & Draine 
(1993) (using the more recent smoothed UV data for silicates) 2 . Here we have chosen the sim- 
plest "standard" values for the input parameters. For example, the models all used the Anders 



2 Available at http://www.astro.princeton.edu/~draine/dust/dust.diel.html 
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& Grevasse (1989) solar abundance set, with the depletion factors for each element given in Ta- 
ble 1. The input ionizing spectrum for the following models was a standard power law continuum 
with (fofa) z^ -1 ' -The effect of varying these and other input parameters will be investigated in 
a following paper (Groves et al. (2001) in preparation). 

To facilitate comparison, and to separate the effects of internal density structure, radiation 
pressure and photoelectric heating by the grains, we have run four families of models with increasing 
ionization parameter Uq (defined at the outer boundary of the ionized slab). In order of increasing 
physical sophistication (and reality) these are: 

1. Isochoric (constant density) dust-free models with heavy element abundances in the gas 
phase which match those of the models including dust. 

2. Isobaric dust-free models. Again these have gas-phase heavy element abundances matching 
the models with dust. These models are isobaric in the sense that the product nT is kept constant, 
but the radiation pressure is not explicitly included. 

3. Isochoric dusty models. 

4. Isobaric dusty models without radiation pressure. These provide a dusty comparison to the 
complete model which shows the difference that radiation pressure makes. 

5. Complete isobaric dusty models, which include the effect of radiation pressure on the density 
structure of the ionized gas. 

To avoid differential effects due to different electron density in the low-excitation cases, in all 
cases the initial density was adjusted so that the density of the region emitting in the [Sn] lines 
was maintained constant at about 550 cm -3 to within 10%. This density was chosen to be typical 
of the electron densities inferred from the [Sn] AA6717/6731 ratio in the NLRs of Seyfert Galaxies. 
The models were terminated at a point where the gas temperature dropped rapidly to below 5000K 
and the emission in the optical forbidden lines is quenched. 

Generally in the literature we encounter only the first or the third of these sets of models. 
The profound difference in the internal structure of the isobaric dusty models including radiation 
pressure compared with a "standard" isochoric model is emphasized in Fig. 2 which shows the 
internal density and temperature structure produced in the dusty, radiation pressure dominated 
slab model with ionization parameter of Uq = 0.5. In this model, the radiation pressure in the EUV 
radiation field is roughly 30 times the initial gas pressure. At the "front" surface of the cloud the 
opacity of the cloud is entirely dominated by dust, and the photoelectric heating by dust is very 
large. As a consequence, the electron temperature in this zone is very high, of order 25000K. The 
electron temperature decreases strongly throughout the model towards the ionization front, and 
this steepens the density gradient in the model. 

This strong density gradient which is set up between the high ionization and low ionization 
regions can provide a physical explanation of the pressure difference between these regions which 
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was required in the A/Mj models of Binette et al. (1996). As expected, the pressure in the low- 
ionization region scales as the radiation pressure at high Uq. Thus, the local ionization parameter 
characterizing the low and intermediate ionization region of the model becomes independent of Uq. 

The emission line ratios produced by the various classes of model can be quite different one 
from another. We have chosen to display here only those which provide a good separation between 
the different classes of model, and so have a diagnostic utility Any of the of line emission ratios 
used by Veilleux & Osterbrock (1987) will give such a separation. 

In figure 3 we plot one of these diagnostic diagrams, the [S n] AA6717 + 673lA/Ha ratio vs. the 
[Olli] A5007A/H/3 ratio. The observational points are from Allen et al. (1999). For the dust-free 
case, both the isochoric nor the isobaric sequence can reproduce the observed data, though the fit 
would be improved with a flatter spectral index. The closest fit is found in the region log Uq ~ —2.8. 
As has been shown by many other authors, only a very restricted range of ionization parameter 
provides a reasonable fit. 

For the dusty models, both the isochoric and isobaric sequences approach the region of the 
observed NLR at high ionization parameter. The dusty, radiation pressure affected isobaric model 
shows very slow variation in the line ratios with U at high ionization parameter and these line 
ratios are in the region of interest. The dusty, isochoric models actually begin to head away from 
the observations as we increase in ionization parameter. 

Figure 4 is a close up of the region of interest, showing only isobaric models. In addition 
to the two isobaric models described above, a dusty, isobaric model without radiation pressure is 
included for comparison. This figure illustrates the self-regulating effect of the radiation pressure 
in producing nearly constant line ratios for Uq > 1CP 2 . While the radiation pressure free case does 
reproduce the observed data quite well, it does so only for a small range of ionization parameter, 
—2.3 < log [To < —1.3. By logfo ~ —1.0 the model has left the region in which the observations 
lie. The complete dusty models, which include the effects of radiation pressure, not only lie in the 
region of the observed data but also display line ratios which stagnate with respect to changes in U, 
over a broad range of values. The implication of this is that, provided Uq S> 10~ 2 , a cloud of given 
density could be moved radially by a large amount without producing any appreciable change 
in its excitation state. Such a behavior is frequently seen, and is interpreted in the "classical" 
photoionization models as indicating a constant ionization parameter. 

The conclusions derived from figure 3 and 4 are reinforced when we replace the [Oni] line 
by the higher-excitation [Nev] line. This is shown in figure 5. Now the range of acceptable 
models is appreciably reduced. From inspection of both Figures 3 and 5 we see that the dusty 
models that include radiation pressure provide acceptable solutions in the range — 1.6 < log Uq < 0. 
The acceptable range of solutions for the dusty isochoric models is much more restricted; —2.0 < 
log Uq < —1.2, and there are no really acceptable solutions for any dust-free model. This illustrates 
a major advantage of the radiation pressure dominated model in that it encompasses a wider 
range of ionization species than any simple isochoric or isobaric model. The reason for this is that 
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as density increases inward, one is effectively spanning a very wide range in the local ionization 
parameter. Furthermore, the extra heating provided by dust photoheating scales with Uq and the 
higher electron temperatures enhance the intensities of the high ionization species in a way that 
the dust free models cannot account for. This effect is further enhanced by dust absorption which 
strongly reduces the luminosity of the hydrogen recombination lines in the high ionization zone. At 
the same time, it enhances the Hen A4686A/H/? ratio to values as high as 0.5-0.6 which is in better 
agreement with observation, and similar to the value obtained in the matter-bounded models of 
Binette et al. (1996). 

If we look only at line ratios which are sensitive to the excitation, then dust-free models might 
seem to provide an acceptable fit. An example is shown in figure 6 which plots the [O I n]/[0 1] line 
ratio against the [Ne v]/[Neiil] ratio. Here similar values of the ionization parameter produce very 
similar values of the [Ne v]/[Ne in] ratio, while the [O I n]/[0 1] line ratio is more strongly affected by 
the temperature in the [O i] - emitting zone. Note that, again the dusty models including radiation 
pressure provide acceptable solutions in the range —1.6 < logf/o < 0. 

The effect of photoelectric heating is best investigated in terms of the electron temperature in 
the [O I n] - emitting region. This is shown in figure 7 which plots the temperature sensitive [O I n] 
A4363A/5007Aratio against [Oiii] A5007A/H/3. For clarity we have omitted the isochoric models, 
since these largely overlap with the isobaric cases. The observed points are drawn from both 
Allen et al. (1999) and Tadhunter, Robinson & Morganti (1989). The dust-free models are clearly 
much cooler than the dusty models and fail to reproduce the observed [O I n] temperatures for any 
ionization parameter. This is the famous "temperature problem" of NLR, and which was discussed 
in detail by Tadhunter, Robinson & Morganti (1989). Photoelectric heating by dust, proposed as 
the solution to a similar problem occurring in planetary nebulae by Dopita & Sutherland (2000) 
seems to provide enough additional heating in NLR to bring the predicted electron temperatures 
into the observed range. Once again, for dusty isobaric radiation-pressure dominated models, 
acceptable models for the [Oiii] A4363A/5007A ratio lie in the range —1.6 < logC/o < 0. 

A final major effect of the addition of dust in these models is to limit the H/3 surface brightness 
as a consequence of the absorption of ionizing photons by dust. Since dust absorption strongly 
dominates over the gas absorption in the regime of high ionization parameter, this limits the 
ionized column (Stromgren column) to a maximum value. For a fixed density or pressure in the 
region emitting the [S II] lines, the surface brightness of the ionized zone in H/? saturates at a 
maximal value at high Uq. By contrast, at low Uq, dust absorption is negligible compared with gas 
absorption, and the H/3 flux scales with the ionization parameter. This behavior is shown in figure 
8. 

If the [S II] density is lowered, then, at at a given Uq the H/3 flux simply scales with the density 
in the low-C/o regime. However, since the dust and the gas absorption terms have a different density 
dependence, at a given Uq in the high-£/o regime, a lowering of the density in the [S II] zone leads to 
a greater fraction of the incident radiation field being absorbed by dust. The result of this is that 
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the saturation value of the H/3 surface brightness is reduced by a greater factor than simply the 
ratio of the densities. The roll-off between the low-t/o and high-C/o, regimes must therefore occur 
at a slightly lower value of Uq than in the high density case. 

The plane-parallel models presented here strictly apply only along the (stagnation) line which 
lies directly towards the ionizing source on the front face of a photoevaporating cloud. Towards 
the sides of the cloud, lateral motions will be developed in the ionized plasma by the transverse 
gas pressure gradient. Around the cloud edges, the radiation pressure force acts radially outwards 
from the central engine. According to equation (2), this drives an accelerating flow radially away 
from the central AGN. Such a radiatively-driven outflow was considered by Binette (1998) in the 
context of the coronal gas. In the case of NGC 1068, new HST STIS data presented by Cecil et al. 
(2001) shows [O in] tails of clouds about 100 pc out from the nucleus being accelerated up to 2000 
km s' 1 or greater, over a distance of about 0.15 arc sec (10 pc). This corresponds to an acceleration 
of order 6 x 10 _4 cm s -2 . Let us now examine whether this acceleration could be due to radiation 
pressure. 

According to the observations of Kraemer & Crenshaw (2000), the density given by the [Sll] 
lines in NGC 1068 is about 2500 cm' 3 . In the high ionization parameter models, this density is 
determined entirely by the radiative pressure, and such a pressure is obtained from a central source 
giving a surface luminosity of ~ 1500 erg cm 2 s _1 at the outer boundary of the photoablating 
cloud. This figure is in excellent agreement with that inferred from the bolometric FIR luminosity 
(Dopita et al. (1998)) ~ 1.5 x 10 11 L . This gives a flux at a distance of 70 pc from the nucleus 
(typical for these clouds) of about 1700 erg cm 2 s -1 . 

If these radiative forces go into accelerating the gas in the ionized tails of the NGC1068 clouds, 
then from equation (4) we would require a flux of order 1.5 x 10 4 erg cm -2 s _1 to provide the 
observed acceleration. The inferred flux is a factor ten lower than this, and this would produce a 
final outflow velocity of order 700km s _1 . The only way around this problem would be to increase 
the force per unit mass of gas by increasing opacity of the dust. This could be done through grain 
shattering, which according to Jones, Tielens, & Hollenbach (1996) will occur whenever grain-grain 
collision velocities are greater than the sound speed in the grain material, about 2 km s _1 . This 
could be achieved with either subsonic turbulence in the gas, or else subsonic radiatively-driven 
grain drift velocities. 

A potential mechanism for grain shattering is as follows. Close to the ionization front, and 
within it, the net grain charge is negative, being dominated by collisions. When the grains are 
charged their motion is strongly coupled to the gas via Coulomb interactions. As the grains are 
advected with the gas flow from the ionization front, the local radiation field becomes stronger 
and the charge on the grains becomes smaller as photoelectric charging becomes more important. 
Eventually, the net charge on an individual grain becomes zero as the collisional and photoelectric 
charging exactly balance. At this point the grain can decouple from the gas flow and radiative 
pressure forces push it back towards the ionization front. The terminal drift velocity that the 
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grain can reach is determined by the gas density, p, and the radiation field intensity, F; Vdriit = 
(F/cp) 1 / 2 . For the parameters computed for NGC 1068 this could amount to ~ 200 kms^ 1 . 
However, long before this velocity is reached, the grain will become charged once more, and the 
grain motion is again locked to the gas motion. Thus, when their charge becomes zero, the grains 
will tend to oscillate about the position where this happens until another (charged) grain strikes 
and shatters them. The shattered fragments have their equilibrium point further out in the ionized 
flow, and of course they are much more vulnerable to the other grain destruction mechanisms we 
have mentioned. In our models for NGC 1068 grains are completely shattered in the region which 
produces [Nev] and [Fevn]. 

The idea that the ionized tails of the cloudlets in NGC 1068 are driven by radiation pres- 
sure acting on the shattered dust seems a reasonable one, but it is clear that the full solution of 
this problem will require a hydrodynamic model which explicitly includes the grain physics and 
destruction processes. We intend to present such a model in a future paper. 

If dust is indeed destroyed at some point in the outflow, the gas will continue to flow out- 
ward from the cloud, but radiation pressure effects become relatively ineffective in determining the 
dynamics. In this case the cloud will become surrounded by a highly ionized coronal region with 
an equilibrium electron temperature of order 8 x 10 4 K. On the basis of continuity of pressure, we 
estimate that in the case of NGC1068, the density of this region is 10 — 20 cm -3 and therefore that 
the effective ionization parameter is of order unity, or somewhat greater. This would be sufficient 
to explain the extreme degree of ionization observed in this coronal region (Kraemer & Crenshaw 
(2000)). 

Finally, let us examine if the hypothesis that dust can be accelerated in the cometary tail of 
the cloud is also consistent with the hypothesis that dust can also be destroyed in entering the 
coronal region. From the acceleration in the tail (6 x 10~ 4 cm s~ 2 ) and the length of the tail (lOpc) 
we can estimate the transit time of material through the tail; 10 4 yr. Along the stagnation line 
in the direction of the source, material flows off the cloud at about the sound speed, 20 km s _1 . 
However, the photoionization model for NGC 1068 shows that the thickness of the ionized zone is 
about 10 18 cm. Therefore, the dusty plasma flows through this region in about 1.6 x 10 4 yr. These 
timescale are quite comparable, which strongly suggests that the length of the [Om] tails and 
the production of the coronal region in NGC 1068 are both moderated by the dust destruction, 
occurring on a timescale of order 10 4 yr. 

4. Conclusions 

In this paper we propose that the NLR and ENLR of Seyfert galaxies can be self-consistently 
modelled by radiation-pressure dominated dusty photoionized regions surrounding photoablating 
dense clouds. By using the code MAPPINGS III, we have shown that these display a strong density 
gradient increasing towards the ionization front, are characterized by strong photoelectric heating 
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by the dust and in which the dominant absorption of the ionizing continuum is due to the dust. 
In these models, the gas pressure close to the ionization front is determined by the radiation field. 
The models provide: 

• A line spectrum consistent with the observations, and which remains approximately invariant 
against changes in Uq that span a range larger than 2 orders of magnitude. 

• A higher electron temperature than can be achieved in standard photoionization models, and 
which could help explain the long standing "temperature problem" of NLR and ENLR. 

• An H/3 surface brightness which reaches a maximum value at high values of the initial ion- 
ization parameter {Jo- 
in addition, they may also offer an explanation for: 

• The co-existence of the low-ionization zones with a extremely high ionization parameter coro- 
nal region. 

• Strong radial radiative acceleration of dusty gas which may explain the broad and blue-shifted 
[O III] profiles seen in many Seyfert galaxies. 

We will systematically explore the properties of such models in future papers. 
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Table 1. Solar metallicity (Z@) and depletion factors (D) adopted for each element. 



Element log(Z ) log(D) 



H 


0.00 


0.00 


He 


-1.01 


0.00 


C 


-3.44 


-0.30 


N 


-3.95 


-0.22 





-3.07 


-0.22 


Ne 


-3.91 


0.00 


Mg 


-4.42 


-0.70 


Si 


-4.45 


-1.00 


S 


-4.79 


0.00 


Ar 


-5.44 


0.00 


Ca 


-5.64 


-2.52 


Fe 


-4.33 


-2.00 
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Fig. 1. — Depiction of our theoretical model showing the photoevaporating, dusty molecular cloud 
(black circle) at the center surrounded by the dusty, photoevaporated plasma, from which both the 
low- (close to the cloud) and intermediate- ionization (close to the coronal region) lines are emitted. 
This plasma is itself surrounded by the dust-free coronal region. The stagnation line indicates the 
region modelled by this paper. Each region (excluding the molecular cloud) is marked with some 
of the strong lines emitted from this part of the structure. The parameters marked by the dashed 
lines indicate typical values for each region for a NLR cloud found in NGC 1068. 



-17- 




Fig. 2. — Temperature and density structure in a radiation-pressure dominated, dusty model with 
U = 0.5 
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Log( [SII]U6725/HcO 

Fig. 3. — Line diagnostic diagram of [O m]A5007/H/3 vs. [S ll]AA6717, 31/Ha showing the output 
from four families of models. These are split into two dust free models (empty markers) and two 
dusty models (crossed markers), with an isochoric model (dashed lines) and an isobaric model (solid 
lines) in each case. The isobaric, dusty model also includes the effect of radiation pressure. Along 
each model line marks indicating log?7o are shown. The solid circles are observations from Allen et 
al. (1999). Note that the dust free cases do not match the observations whereas both dusty models 
provide reasonable solutions. 
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Log( [SII]AAtf725/Ha ) 

Fig. 4. — A close-up of the region of interest in figure 3 showing only the isobaric models. As well 
as the dust free isobaric model (empty circles) and the dusty, radiation pressure affected model 
(crossed circles) from the previous figure, a dusty model without the effects of radiation pressure is 
included (dashed line with dotted circles) . Note the model without P ra d does match the observations 
well, but the line ratios do not stagnate with changes in U like the complete dusty model. 
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Fig. 5. — Line diagnostic diagram of [Ne v]A3426/H/3 vs. [S ll]AA6717, 31/Ha with the same symbols 
as figure 3. Note that in consideration with 3 that only the dusty models provide reasonable 
solutions the isochoric case with —2.0 < log Uq < —1.2 and the isobaric case with —1.6 < log L^o < 0. 
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Log( [O m]5007/[O I]6300 ) 

Fig. 6.— Line diagnostic diagram of [Ne v]A3426/[Ne in] A3968 vs. [0 in] A5007/[O i] A6300, with the 
same symbols as figure 3. Note that these line ratios are sensitive to excitation, and with inspection 
of this and the previous two figure it can be seen that again the isobaric dusty model provides a 
reasonable solution between —1.6 < log Uq < 0. 
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Fig. 7. — Line diagnostic diagram of [O m]A5007/H/3 vs. [O lli]A4363/A5007 with two cases: a dust 
free isobaric model (empty circles) and a dusty, isobaric with radiation pressure model (crossed 
circles). The marks on the model lines indicate logt/Q- The solid circles are observations drawn 
from both Allen et al. (1999) and Tadhunter, Robinson k Morganti (1989). Note that with the 
consideration of figures 3-7 that the isobaric, dusty, radiation pressure affected models provide the 
most acceptable solutions with —1.6 < log[/o < 0. 
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Fig. 8. — H/3 flux vs. ionization parameter, showing two groups of models. The two groups shown 
differ in density, with the crosses indicating a density at the back of the cloud of ~ 550 cm -3 and 
the circles, ~ 55 cm~ 3 (as labelled). Within each group the thin lines show a 'standard' (isochoric) 
model without dust, and the thick lines the result from our radiation-pressure dominated, dusty 
model, showing the limiting of the H/3 flux. 



